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We report on the analysis of the root mean square curvature as a function of the numerical resolution for a 
single reconnection of two quantized vortex rings in superfluid helium. We find a similar scaling relation as 
reported in the case of decaying thermal counterflow simulations by L. Kondaurova et al. There the scaling 
was related to the existence of a Kelvin-wave cascade which was suggested to support the L’vov-Nazarenko 
spectrum. Here we provide an alternative explanation that does not involve the Kelvin-wave cascade but is due 
to the sharp cusp generated by a reconnection event in a situation where the maximum curvature is limited by 
the computational resolution. We also suggest a method for identifying the Kelvin spectrum based on the decay 
of the rms curvature by mutual friction. Our vortex filament simulation calculations show that the spectrum of 
Kelvin waves after the reconnection is not simply n{k) == with constant T]. At large scales the spectrum 
seems to be close to the Vinen prediction with TJ = 3 but becomes steeper at smaller scales. 


I. INTRODUCTION 

The decay of quantum turbulence in the limit of zero tem¬ 
perature is one of the key topics that seeks an answer. The 
Kelvin-wave cascade is one promising way of transferring en¬ 
ergy from intervortex distances to smaller scales where the 
classical Kolmogorov cascade is impossible There the 
weak nonlinear coupling between different scales of heli¬ 
cal distortions (=Kelvin waves) transfers the energy to ever 
smaller scales until dissipation may become possible, e.g., via 
phonon emission. 

Reconnections and interactions between neighboring vor¬ 
tices are expected to work as a drive that transfers energy from 
large scales and feeds the Kelvin-wave cascade. At finite tem¬ 
peratures the weak Kelvin cascade is easily suppressed by the 
mutual friction dissipation^’®. In order to understand the de¬ 
cay processes in the zero temperature limit at scales around 
and smaller than the intervortex distance, we must be able to 
distinguish between the roles of the Kelvin-wave cascade, the 
Kelvin waves directly induced by the reconnection events, and 
the mutual friction dissipation. 

Here we concentrate on the Kelvin waves generated by a 
single reconnection event and investigate whether a Kelvin 
cascade originates from the reconnection event. We also illus¬ 
trate how Kelvin waves are damped by mutual friction. Our 
analysis at zero temperature reveals that soon after a recon¬ 
nection the vortex becomes filled with Kelvin waves of differ¬ 
ent scale and that the presence of the shortest scales is limited 
only by the resolution. However, in accordance with earlier 
work^, the evidence for the existence of the Kelvin cascade 
remains weak since the observations can be explained by a re¬ 
distribution of Kelvin waves that are initially “packed inside 
the reconnection cusp”. 

Our paper is organized such that in Sec. II we review the 
relation between the rms curvature and the Kelvin spectrum. 
We also take into account the numerical contribution from a 
single sharp cusp. In Sec. Ill we consider a straight vortex 
with Kelvin waves and show how the decay of rms curvature 
by mutual friction can be used to identify the Kelvin spectrum. 
We additionally illustrate that the mean curvature can strongly 
increase without changes in the Kelvin spectrum. In Sec. IV 
we consider the actual reconnection event and use the results 


from previous sections to interpret these results. 

II. MODEL 

In superfluid ^He the vorticity appears in the form of quan¬ 
tized vortices, with circulation quantum K = hjnn k, 0.0997 
mm^/s, whose core size is of the order of 1 A. Ignoring com¬ 
pressibility effects, the vortex dynamics can be modeled by 
using the vortex filament model where the vortices are con¬ 
sidered as line defects and the superfluid velocity field is 
given by the Biot-Savart law. The filament model should 
be a good model for ^He-II where the core radius is five to 
six orders of magnitude below the typical experimental inter¬ 
vortex distance. The model is extensively described, e.g., in 
Ref. [8]. Numerically the vortex, a three-dimensional curve 
s(^,f), where ^ is the length along the vortex and t is time, is 
described by a sequence of points whose motion is here fol¬ 
lowed by the fourth order Runge-Kutta method. Our numerics 
is described in more detail in Ref. [7], where we analyzed a 
similar situation of reconnecting vortex rings but concentrated 
on the dissipation due to mutual friction. 

A. Kelvin spectrum and rms curvature 

The identification of helical Kelvin waves in a vortex tan¬ 
gle is a delicate problem®. A proper identification is currently 
possible only in a few rare cases. In the case of a straight vor¬ 
tex, taken to be along the z axis, and when the perturbations 
can be presented as w(z) = x{z) -f iy(z), the Fourier transfor¬ 
mation of w{z) determines the Kelvin-wave amplitudes w^. 
The Kelvin occupation spectrum is typically defined as 

n{k) = \wk\^+\w^k\^ ,k>0. (1) 

One should note that here the k vector is a one-dimensional 
vector that defines the wavelength for the Kelvin waves and 
should not be confused with the amplitude of the three- 
dimensional k vector that is typically used when writing the 
energy (velocity) spectrum. 

To provide a rough measure of the Kelvin waves the dis¬ 
tribution of curvature radii can be extracted from a vortex 
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filament calculation. However, in the general case, to de¬ 
rive the Kelvin spectrum from the curvature spectrum leads 
to an insolvable inverse problem. In case of a straight vortex, 
and in the limit of small Kelvin amplitudes, these two spec¬ 
tra are closely related. Unfortunately, the curvature spectrum 
is typically much too noisy and no proper spectrum can be 
determined from the curvature data. Therefore, the analysis 
in Ref. [6] was restricted to the calculation of the root mean 
square (rms) curvature, Cnn.,: 

= ( 2 ) 

where c(^) = |s"(^)| is the local curvature at and L is 
the total vortex length. For small Kelvin amplitudes (more 
precisely when |w'(z)| = \dw/dz\ <C 1) and when the en¬ 
ergy is taken to be proportional to the vortex length [local¬ 
ized induction approximation, (LIA)], the energy spectrum 
E{k) oc k^n[k) for Kelvin waves with wave number k can be 
related to the rms curvature by® 

k^Eik)dk. (3) 
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scales smaller than the intervortex distance transports energy 
towards the core scales where it can be dissipated. For this 
to be true, at sufficiently low temperatures the Kelvin-wave 
cascade causes all our numerically traceable scales to be oc¬ 
cupied by Kelvin waves, and the rms curvature becomes deter¬ 
mined by the numerical resolution. By repeating simulations 
with different resolutions one should be able to distinguish 
between different theories. 

This is precisely what was done in Ref. [6] for decaying 
thermal counterflow, where the scaling for Crms weakly sup¬ 
ported the explanation in terms of the Kelvin-wave cascade 
with L’vov-Nazarenko spectrum. The complicated situation 
was exemplified by the fact that the resolution could only be 
changed by a factor of 2.7, and no firm proof was obtained. 

In the following we consider a simpler case; a single recon¬ 
nection of two initially linked vortices. In this case we are able 
to change the resolution by a greater amount. However, first 
we consider how the rms curvature is affected when a vortex 
configuration contains a sharp cusp where the curvature may 
diverge (or at least yield values close to the inverse core size). 

B. Sharp cusp and its effect on rms curvature 


The term A = ln(£/ao) depends only weakly on the intervor¬ 
tex distance i and on the core radius aq. Here the Kelvin spec¬ 
trum is assumed to be valid from up to some cutoff scale 
^max, above which the spectrum is suppressed (either due to 
numerical cutoff or, e.g., by suppression due to mutual fric¬ 
tion). 

Since the predictions for the Kelvin-wave spectrum vary 
from E{k) l/k, in the case of strong wave turbulence 
with the Vinen spectrum^**, to E{k) oc with the L’vov- 
Nazarenko spectrum^, the rms curvature is not much affected 
by the scales near kmin. The dominant contribution originates 
from scales near kmax and may therefore be limited by the nu¬ 
merical resolution of the calculation. 

The different proposals for Kelvin-wave spectra by Vinen^® 
(V), Kozik-Svistunov^ (KS), and L’vov-Nazarenko^ (LN) in 
terms of Crms are asymptotically given by® 

^mis “ *^^max 

= <I>(£k^ax)^/^ (4) 

Generally, when the Kelvin spectrum n{k) k^'^, the rms cur¬ 
vature is given by Crms The intervortex spacing £ 

limits the smallest possible k values, and the term <I>, given by 

! - rkmsa. 

^ = ^AnE/KK^, E= E{k)dk, (5) 


As soon as one introduces a reconnection in the filament 
model one creates a sharp cusp in the vortex configuration. 
A simplified analytical model for a cusp where the curve 
has a simple discontinuity in the derivate, being smooth else¬ 
where, does not necessarily change the rms curvature by a 
large amount. However, numerically the effect can be much 
more dramatic. Without additional smoothing, the maximum 
curvature at the cusp is the inverse resolution, and the region 
for this curvature peak is of the order of the point separation 
(A^res)- Smoothing does not necessarily solve the problem, 
because it only widens the region and restricts the maximum 
value. 

If we split the integral into the “smooth part”, far from the 
cusp, and into the region that is strongly affected by the nu¬ 
merics, one obtains that 
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Here the numerical factor b is of the order of unity. It should 
depend weakly on the numerical details and, e.g., on the initial 
angle between the vortices before the reconnection. In the 
limit of high resolution the last term becomes dominant, and 
the rms curvature due to numerics becomes 


c 


N 

rms 



k 


1/2 
max 5 


(7) 


takes into account the fraction from the total energy E (per 
unit length, per unit mass) related to Kelvin waves. Because 
some numerical factors of the order of unity have been ne¬ 
glected in the above derivation, one should consider <I> only as 
an estimate. 

We expect that at low enough temperatures the Kelvin cas¬ 
cade becomes the dominant dissipation mechanism, which at 


where kmax = 2;r/A^re.s is determined by the resolution. 

If we consider the analytical model (which is based on ear¬ 
lier numerical work) for the vortex shape before the reconnec¬ 
tion [11], we get further support for Eq. (7). This estimate 
for the prereconnection vortex shape indicates that the cusp 
sharpens until the minimum distance reaches the core size (af¬ 
ter which the filament model is not valid any more) such that 
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eventually the region, where the curvature is of the order of 
the inverse core radius, is also of the order of the core diame¬ 
ter. Simulations with the Gross-Pitaevskii model indicate that 
the curvature does not increase much from this limit imposed 
by the core radius and that it also has a similar magnitude af¬ 
ter the reconnection^^. Therefore, Eq. (7) should originate al¬ 
ready from the prereconnection dynamics where k^ax is given 
by the core size. 

Even if the cusp contribution to Crms grows slower with 
kmax than the cascade predictions in Eqs. (4), it may give the 
dominant contribution when the Kelvin amplitudes are small 
{<P ^ 1) and when the resolution is not high enough, as shown 
below. Its contribution also increases with increasing number 
of reconnections. 

In addition to the rms curvature, we also use the mean cur¬ 
vature Cave in our analysis; 

Cave = (8) 

Based on similar arguments as used for the rms curvature in 
Eq. (6), the contribution from the cusp region should result in 
a constant contribution that does not increase as the resolution 
is improved, i.e., 

Cave = |/ c{^)d^ + j, (9) 

^ J smooth ^ 

where b is another constant of the order of unity. Because 
Cave is sensitive to the phases of the Kelvin waves it cannot be 
directly related to the Kelvin spectrum. 


we obtain by using the small amplitude estimation, where 
c(^) ~ |w"(z)| and ~ dz, that 
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Here the vortex length L = Lit) k, depends on time only 
weakly. By taking kmax —°° and approximating the sum with 
an integral we obtain, for Tj < 5, the asymptotic decay curve 
for the rms curvature as 
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where r(fl) is the gamma function. If for briefness we write 
Crms = D(at)^^, then the decay exponent = (5 — J7)/4 
is simply given by the Kelvin spectrum. Therefore, if the 
Kelvin-wave cascade is negligible during the time window 
when the decay occurs, the above law can be used to identify 
the initial Kelvin spectrum, as shown below. 

If one wants to improve the above approximation, Eq. (12), 
at small times by taking into account the finite resolution, 
one should only replace the gamma function with the in¬ 
complete (lower) gamma function y(a,x) as r((5 — T])/2) —>■ 
7((5 — 77)/2, 2atm^^^/Til). Here Wmax is the largest Kelvin 
mode present, which is limited by the resolution. At large 
times Eq. (12) fails when only a few lowest modes are present 
(dominant). Eventually the decay is simply given by the ex¬ 
ponential decay of modes m = ±1. 


III. KELVIN WAVES ON A STRAIGHT VORTEX 

In order to better understand simulations where Kelvin 
waves are produced by a single reconnection event, we first 
consider two simplified sets of simulations for a straight vor¬ 
tex where the determination of the Kelvin spectrum is less 
challenging. These calculations model the time development 
after the reconnection. In both cases we occupy the straight 
vortex, taken to be along the z axis, with a known Kelvin spec¬ 
trum parametrized as 

n{m) = \w,n\^+ \w-m\^ , m = 1,2,3,... . (10) 

Here A determines the amplitude of the Kelvin spectrum, and 
the Wm’s can be obtained by Eourier transformation of w(z), 
provided that the configuration remains single valued. At the 
instant when the single valuedness breaks, \dw/dz\ diverges. 
The k vector k = 2nmlL^ is given by the mode number m and 
period L^. In all these calculations we have fixed the reactive 
mutual friction parameter a' = 0. 

We can make an analytical approximation for the time de¬ 
velopment of Crms if we assume that the Kelvin-wave cascade 
is absent such that every Kelvin mode decays exponentially as 
Wm{t) = where Tam = Tii/(aOT^). With Tu we 

denote the decay time for mode m = 1 when a — \. We ne¬ 
glect the weak logarithmic m dependence in Tam- By assum¬ 
ing that initially Wm(0) has a spectrum defined by Eq. (10), 


A. Decay of rms curvature via mutual friction 

In our first sample case, we use different Kelvin spectra 
with different amplitudes, and observe how the rms curvature 
decays by the dissipative mutual friction (x{T). We have addi¬ 
tionally assumed that initially the phases of the Kelvin waves 
are random and that |wm| = |w_m|. Eigure 1 summarizes our 
results in case of A = O.OIL, and a =0.1. The number of 
points used is 1024 and L~=\ mm. 

By repeating the calculations with different values of a, we 
observe that the decay curve Crms = D{at)^P is universal, as 
predicted by Eq. (12), and the terms D and p are not essen¬ 
tially affected by the mutual friction parameter a. Also, the 
exponent j3 is almost solely given by the Kelvin spectrum ex¬ 
ponent 7], as shown by the lower inset of Eig. 1, where we 
have used a time window 10^^ s < at < 10^^ s for the fit. 
The lower (red) curve is calculated for a relatively large am¬ 
plitude A/L, = 0.1, where the energy related to Kelvin waves 
is of similar order as the energy related to the straight vortex. 
The upper (blue) curve is appropriate for the small amplitude 
limit A/L; = 0.001. Eurther expansion of the amplitude win¬ 
dow does not essentially change the results. 

The amplitude D is mainly determined by the Kelvin ampli¬ 
tude A. This is illustrated in the upper inset of Eig. 1 using the 
L’vov-Nazarenko spectrum. The small time behavior, which 
originates from the smallest scales, is sensitive to the resolu¬ 
tion and is consistent with the saturation if the full Eq. (11), or 
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FIG. 1. (Color online) Decay of Crms for different initial Kelvin spec¬ 
tra on a straight vortex with periodic boundary conditions and period 
L, = 1 mm. The number of points is 1024 and a = 0.1. In the 
main panel the amplitude is A/Lj = 0.01, see Eq. (10). The different 
spectra are from top to bottom defined by T] = 3.0 (red, Vinen), 3.4 
(magenta, Kozik-Svistunov), 3.666... (black, L’vov-Nazarenko), 4.0 
(blue), and 4.4 (cyan), respectively. The dashed lines are from the 
analytical model, Eq. (12). The upper inset shows the Crms for the 
LN spectrum but for different amplitudes: from top to bottom A/L^ 
= 0.3, 0.1, 0.03, 0.01, 0.003, and 0.001, respectively. The lower in¬ 
set relates the decay exponent j6 in Crms ^ to the exponent T] 
in the Kelvin spectrum of Eq. (10). The upper (blue) curve is for 
small amplitudes with A/Lj = 0.001, and the lower (red) curve is for 
rather large amplitudes with A /L- = 0.1. The thin dashed line is the 
exponent appearing in the analytical result [see Eq. (12)]. 


the approximation described below Eq. (12), is used. There¬ 
fore, the power law can be validated also at early times simply 
by improving the resolution. 

If we use the value Tn Ri 0.22 s from the calculations where 
only the mode m = 1 is present and approximate that L k, L-, 
we obtain nice quantitative agreement with Eq. (12), not only 
for the exponent, but also for the absolute amplitude. This 
is illustrated in Eig. 1 with the dashed lines. The general 
tendency for deviations upwards from the theory prediction, 
Eq. (12), originates from the omitted logarithmic ln(m) depen¬ 
dence in T,na- The deviations downwards at large times (which 
is visible in Eig. 1 for t] = 4.4 and 4.0) can be understood to 
originate from the finite period which results in that eventu¬ 
ally only the longest Kelvin waves with wavelength equal to 
L, persist and decay exponentially. Equation (11) models this 
limit accurately. 

The decay curve for Cnns is not sensitive to the initial phases 
of the Kelvin waves, as expected from the derivation where the 
phases are not present. We have also numerically checked that 
on the scale of Eig. 1 one cannot resolve two different decay 
curves with different initial sets of random phases. 

If we analyze the decay of the different Kelvin modes in 
Eig. 1, we may verify that the different Kelvin modes actually 
rather well exhibit exponential decay where the time scale is 
determined by a and k. There are additionally fluctuations 
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FIG. 2. (Color online) Time development of the “test-cusp” vortex 
at r = 0 that initially hosts the Kozik-Svistunov spectrum. The up¬ 
per panel illustrates the time development of the local curvature at 
times 0 (black), 0.5 (red), 1 (blue), 2 (green), 10 (cyan), and 100 
ms (magenta). The lower left panel shows the corresponding vortex 
configurations, and the lower right panel contains the Kelvin spectra 
where k = iKnijL- is given by the mode number m. 

that are due to interactions between different modes. Also, 
the decay only continues for about four to five orders of mag¬ 
nitude, which is likely due to numerics. We have not checked 
the situation for very small a values because it requires a large 
enough time window. This would be important if one seeks 
the Kelvin cascade, but if the intention is to determine the ini¬ 
tial Kelvin spectrum one is not required to use extremely small 
values of mutual friction. Actually, it is much more conve¬ 
nient to use a large enough a in order to minimize the effect 
from the cascade during the decay of Crms- 

Therefore, if the characteristic curvatures are much larger 
than the inverse of the intervortex distance l/£, the decay 
of Crnis should provide a potential way to identify the Kelvin 
spectrum, perhaps also in more complicated tangles. 


B. Cusp built by Kelvin waves 

In our second example case we illustrate how a cusp, where 
the curvature is strongly peaked, can be built on a straight 
vortex using a Kelvin spectrum defined by Eq. (10). We then 
follow its dynamics both at T = 0 and T >0. 

The phases of the Kelvin waves are typically considered to 
be random, which results in a wiggly looking vortex. How¬ 
ever, by specially organizing the phases one may generate a 
vortex that is otherwise smooth but contains a single cusp (per 
period). In Eig. 2 we have built a “test-cusp” by setting the 
same constant phase for all the Kelvin waves using the Kozik- 
Svistunov spectrum. The period is L, = 1 mm, and the num¬ 
ber of points is 1024. All the Kelvin waves have rather small 
amplitude A/L^ = 0.01, but since the phases are the same, 
|w'{z)| Ri 0.5 near the cusp. Therefore, the theory predictions 
made in Sec. II are not necessarily valid, e.g., the rms curva- 
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FIG. 3. (Color online) Decay of the rms curvature when the ini¬ 
tial configuration is a specially arranged sharp cusp with Kozik- 
Svistunov spectrum, see Fig. 2. The different curves correspond to 
different values of mutual friction: a = 0.0001 (blue), 0.001 (cyan), 
0.01 (green), 0.1 (magenta), and 1 (red). The inset shows the time 
development for Crms, and in addition at a = 0 (black) both for Cms 
(solid) and Cave (dashed), plotted on a linear time scale. 

ture grows slower than the Cnns ^max prediction. 

The upper panel of Fig. 2 shows the time development of 
the local curvature at T = 0. It seems that with increasing 
time more small scale structure is developing because Kelvin 
waves spread from the cusp and eventually fill the whole vor¬ 
tex. However, if one investigates the corresponding Kelvin 
spectra, shown in the lower right panel of Fig. 2, the amount 
of Kelvin waves n(k)dk) is not essentially changed be¬ 

cause the spectrum is still the same (initial) Kozik-Svistunov 
spectrum. Deformations occur only near the resolution limit, 
where the numerical errors are largest. 

As the cusp relaxes, the phases of the Kelvin waves ran¬ 
domize because all the Kelvin waves have different propaga¬ 
tion velocity along, and different rotation velocity around, the 
z axis. This results in that the mean curvature, which initially 
is much smaller than the rms curvature, grows and soon has a 
similar value as Cnns- This is illustrated in the inset of Fig. 3, 
where we have analyzed the decay of the rms curvature by mu¬ 
tual friction, for different values of a. The main panel shows 
that the decay follows the universal power law with exponent 
j3 = 0.4, derived above for the KS spectrum (Fig. 1). 


IV. RECONNECTION OF TWO LINKED VORTEX RINGS 
A. Numerical setup 

We consider two vortex rings with radius R = \ mm, ini¬ 
tially linked and oriented perpendicular to each other. The 
initial distance between the two vortices is O.ll?, causing that 
at r = 0 they will reconnect approximately at 1.696 s. For 
dynamics, see Ref. [7]. Our reconnection method is rather 
standard: We reconnect the vortex segments when the mini¬ 


mum distance becomes smaller than O.SA^min. where Ai^min 
is the minimum point separation tolerated. The maximum 
point separation allowed is A^max = 2A^niin- Additionally, we 
require that the vortex length must decrease, which approxi¬ 
mately means that the reconnection event is dissipative. 

The time resolution is chosen such that our numerical time 
step is kept smaller than the time scale related to the small¬ 
est resolvable Kelvin waves. This implies that our time step 
At oc (A^jnin)^> resulting in that the better the spatial resolu¬ 
tion, the more difficult it becomes to cover the same overall 
time window as in the lowest resolution run. Another numeri¬ 
cal challenge is caused by the fact that the computational work 
per one time step grows as N'^ (1/A^niin)^. where N is the 

number of points used to describe the vortex. 


B. Finite mutual friction 

Before concentrating on the zero temperature limit we show 
how finite mutual friction affects the vortex dynamics and 
how it damps the Kelvin waves generated by the reconnec¬ 
tion cusp. In Fig. 4 we have plotted the rms curvature for 
one particular resolution, defined by Ai^min = 0.0025 mm, for 
seven different values of the mutual friction parameter a (in 
all cases a' = 0). 

For a ~ 1 the instant of reconnection depends essentially 
on a because the vortex rings shrink more and therefore they 
also move faster. The relaxation rate for Cnn.s is strongly af¬ 
fected by the mutual friction. However, one may notice that 
the peak value, just after the reconnection, is not substantially 
affected by a. The changes are mainly due to the numeri¬ 
cal algorithm used for the reconnection. This causes that the 
distance between the neighboring vortices is not exactly the 
same when the reconnection is made, and therefore the sharp¬ 
ness/angle of the cusp may slightly vary. 

If we compare these results with the model system, 
Sec. Ill A, we may conclude that the results are consistent with 
the idea that the amplitude of the Kelvin wave decays expo¬ 
nentially with time scale T l/{ak^). These Kelvin waves 
are generated by the reconnection, and the a scaling shown 
in the main panel of Fig. 4 is consistent with the idea that the 
Kelvin-wave cascade is absent. If, at low enough tempera¬ 
tures, the Kelvin cascade would become important (compared 
with mutual friction), then the same scaling would not work 
anymore because the cascade should continuously generate 
new small scale Kelvin waves and Cms should decay slower 
(or perhaps even grow). Note that a similar scaling relation is 
obtained and plotted for the mutual friction dissipation power 
in Fig. 5 of Ref. [7]. 

The absence of the Kelvin cascade is additionally supported 
by the decay curve, plotted with the black dash-dotted line in 
Fig. 4. Here the decay is calculated using a = 0.01, but the 
initial configuration is taken from the zero temperature calcu¬ 
lations with a = 0 one second after the reconnection. Since 
the decay follows exactly the same curve as the standard a 
= 0.01 case, the two spectra are likely the same. Therefore, 
even during that one second period, which corresponds to the 
Kelvin period with a wavelength of the order of R = 1 mm, the 
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FIG. 4. (Color online) Time development of the rms curvature for 
different values of the mutual friction parameter a (and a' = 0). In 
the main figure time is measured from the reconnection instant and 
scaled by a, and on the vertical axis the initial curvature Crjns(O) = 
1/1? is subtracted from the data. Different colors are used for a = 1 
(red), 0.1 (magenta), 0.01 (green), 0.001 (cyan), 0.0001 (blue), and 
0.00001 (light green), respectively. The corresponding maximum 
times simulated are 1, 9, 10, 12, 20, and 30 s, respectively. The black 
dash-dotted curve is for a = 0.01 where the initial configuration is 
taken from T = 0 simulations, 1 s after the reconnection, and the time 
(which covers a time window of 5 s) is measured from this instant. 
The curves are calculated with a resolution of AiJmin = 0.0025 mm, 
while for a = 0.01 results with = 0.00125 mm (dashed, green) 
and 0.000625 mm (dash-dotted, green) are also shown. The power- 
law behavior is illustrated by plotting the Crms (ar)^^ asymptotics 
with P = 0.5 and 0.38 (dashed) that correspond to T] = 3 and 3.48, 
respectively. In the lower inset linear scales are used without any 
scaling. The color coding is the same as in the main figure. The 
case a = 0 is additionally shown by the black curve. The upper inset 
is for a = 0.01 with five different resolutions indicated on the plot, 
together with the j6 =0.38 asymptotic (dashed line). 


interactions between different Kelvin modes have not altered 
the spectrum which was directly generated by the reconnec¬ 
tion at r = 0. 

The observed decay exponent in our reconnection calcula¬ 
tion (see Fig. 4) at late times (large scales) is close to j3 = 0.5, 
which corresponds to the Vinen spectrum (rj = 3) for the 
Kelvin waves. This has been suggested, e.g., in Ref. [10]. At 
early times, which correspond to small scales, the decay slope 
is slower, and one obtains a reasonable fit using j3 = 0.38. Us¬ 
ing the analytical model in Eq. (12), this corresponds to rj = 
3.48. In the upper inset the decay behavior with a = 0.01 
is illustrated using five different resolutions. The two high¬ 
est resolutions are also included in the main panel of Fig. 4. 
The slope at early times is close to the predictions by Kozik- 
Svistunov^ with Tj = 3.4, but here the spectrum is generated 
directly by the reconnection, without the need for a cascade. 

The difficulty in finding a good fit with fixed j5 , that would 
cover the whole time window for the decay of Crms, may in¬ 
dicate that the Kelvin spectrum is not necessarily given by 
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FIG. 5. (Color online) Time development of the rms curvature at 
T = 0 for five different resolutions where the point separation is kept 
between Ai^nijn... 2Ai^jj,in. For the two highest resolution runs, times 
well before the reconnection were omitted and the calculations were 
continued using the configuration from the run with one step lower 
resolution where the configuration was still smooth. The left inset 
shows the detailed behavior near the reconnection event when Ai^[j,in 
= 0.000625 mm. The right inset shows the scaling relation for the 
root mean square curvature after the reconnection as a function of 
resolution plus the fit Cims (solid line) and the predictions 

from the L’vov-Nazarenko spectrum (lower dashed) and from the sin¬ 
gle cusp model (upper dashed). 


Eq. (10) with constant Tj. It seems that the smaller the scale 
the bigger is tj. This is consistent with the idea suggested by 
Nazarenko in Ref. [13], where he argued that the reconnection 
drives not only length scales near the intervortex distance but 
also smaller scales. Therefore, the spectrum excited by the 
reconnection would depend on scales. Now these simulations 
illustrate that the spectrum is not simply n{k) oc that can 
be associated with a sharp bend. 


C. Zero temperature results 

In the zero temperature limit, where the mutual friction de¬ 
cay time for the single Kelvin wave diverges, Cnn.s does not 
essentially decay after the reconnection kink has appeared but 
remains sensitive to the numerical resolution. Figure 5 illus¬ 
trates the time development of the calculated rms curvature 
when the resolution is changed by a factor of 16 while a = 0. 
The inset illustrates that the curvature starts to increase al¬ 
ready before the reconnection due to nonlocal interactions. 
This tends to strongly deform the vortex shape near the re¬ 
connection point. The noticeable jump at the instant of recon¬ 
nection appears due to the appearance of the sharp cusp which 
is sharpest immediately after the numerical reconnection pro¬ 
cess. Therefore, very locally the curvatures reach the values 
determined by the resolution. 

After the reconnection the rms curvature may slightly in¬ 
crease and later decrease, but very soon after the reconnection 
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FIG. 6. (Color online) Time development of the mean curvature 
for the same five different resolutions which were used to plot the 
rms curvature in Fig. 5. The inset shows a similar scaling relation 
as that for Crms- The solid line is a fit, Crms 

points, and the dashed lines are again the predictions from the L’vov- 
Nazarenko spectrum (lower dashed) and from the single cusp model 
(upper dashed). 


Cnns remains constant. This is consistent with the LIA approx¬ 
imation where Crms is one of the several constants of motion 
and where cascade formation is not allowed. We have con¬ 
tinued our Ai^niin = 0.005 mm calculation up to 100 s, and in 
addition to the small ~1% decay of Cnns during this time; it 
possesses oscillations with amplitude ~1% and period ~ 20 
s. This makes it difficult to say accurately from shorter simu¬ 
lation runs whether the curvature actually decays, or whether 
the initial small decay is just part of these oscillations. 

If we plot the approximately constant value of Cnns after the 
reconnection as a function of the resolution one may observe 
a scaling relation, similar to the predictions above. This is 
shown in the inset of Fig. 5. The fitted exponent is 0.56, which 
is closer to the prediction for a simple cusp than the one orig¬ 
inating from the Kelvin cascade with L’vov-Nazarenko spec¬ 
trum. If we want to relate the exponent to a particular spec¬ 
trum, the value corresponds to rj = 3.88. Within the cusp 
model, if we use the minimum point separation instead of 
A^res in Eq. (7), we obtain bK.1. This is consistent with our 
prediction (note that our configuration has two cusps). 

In Fig. 6 we have plotted the time development of the mean 
curvature. One may observe that the large sudden jump, 
present for the rms curvature, is now missing at the instant 
of reconnection. Rather, the mean curvature starts to increase 
after the reconnection and finally stabilizes to a value very 
similar to Crms- Therefore, also the Cave satisfies a similar scal¬ 
ing relation as Crms- 

The time scale for reaching the final value in Fig. 6 is deter¬ 
mined by the resolution, and the increase in Cave appears be¬ 
cause small scale Kelvin waves travel from the reconnection 
site and quickly fill the whole vortex. The better the resolu¬ 
tion, the faster the smallest Kelvin waves travel, and therefore 
the time scale for reaching the “steady-state” value for Cave is 




FIG. 7. (Color online) Time development of the local curvature at 
T = 0 along the vortex after the reconnection of two rings. The 
reconnection results in two strong peaks at locations ~ 0 and 
^ ~ 0.5L. The upper panel is for the resolution - 0.00125 mm, 
and the lower panel is for twice better resolution, = 0.000625 
mm. The black curve is immediately after the reconnection, and 
other curves are 1 ms (red), 10 ms (blue), and 20 ms (green) af¬ 
ter the reconnection. Note the logarithmic vertical axis and that the 
maximum speed of the spreading Kelvin waves is limited/defined by 
the resolution. 


directly proportional to the point separation used. This redis¬ 
tribution of the Kelvin waves is illustrated in Fig. 7, where the 
local curvatures are plotted at different times. 

One might argue that the increase in Cave proves the exis¬ 
tence of the Kelvin cascade because it seems that more small 
scale structure appears. This is possible, but it is not the only 
possible explanation. First, scaling relations similar to those 
in Eq. (4) cannot be easily derived for Cave from the Kelvin- 
wave spectrum. Second, the mean curvature can be strongly 
increased by simply adjusting the phases of the Kelvin waves 
without changing the Kelvin spectrum, i.e., without changing 
the amplitudes of small scale Kelvin waves. This is illustrated 
by our “test-cusp” in Sec. IIIB, see Fig. 3. In the small am¬ 
plitude limit, which is not necessarily a proper assumption for 
a cusp, the phases of Kelvin waves do not change Crms, as 
shown by Eq. (3). Instead, the increase in the average curva¬ 
ture Cave in Fig. 6 after the reconnection appears as the phases 
“randomize” from their initial values at the instant of recon¬ 
nection. 

The above phase randomization can be understood by con¬ 
sidering a vortex with a shape of a triangle wave whose mean 
curvature is zero. It has a Fourier representation with only odd 
indices k = 1,3,5,... whose amplitudes go like wt 1/k^, 
corresponding to the Kelvin spectrum E[k) oc \j]^ (equal to 
the Kelvin occupation spectrum n]^’^ 

from the FN spectrum. The Kelvin waves present in this kind 
of configuration have a different frequency around the “av¬ 
erage vortex direction” (the direction which is realized when 
the amplitude of the triangle wave goes to zero) and differ¬ 
ent propagation velocities. Therefore, the vortex has soon a 
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rather different shape with a mean curvature much larger than 
initially. This does not require any cascade. 

One may still wonder why the Kelvin spectrum, determined 
from the decay curve at small scales where rj ^ 3.5, does 
not show up in the scaling relation for Crms, which is close to 
the cusp model (or alternatively would correspond to Kelvin 
spectrum with 77 « 3.9). One possible reason for this differ¬ 
ence might be that the amplitude of the Kelvin spectrum is 
too small, and therefore the cusp contribution still dominates 
when calculating the rms curvature. If we compare Eqs. (4) 
and (7) we find, by setting ^max = that these esti¬ 

mates give the same value for Cnns when 

(13) 

where 7 = {5 — rj)/2 is the general exponent in Eq. (4) for 
^max- We may crudely approximate that the energy associated 
to Kelvin waves is related to the increase in vortex length, 
which is around 1 % in our example case^^, and therefore 
<I> ft! 0.1. If we use b = 1 from our fit, L ft AnR, and addi¬ 
tionally approximate that I = R = I mm, we obtain that the 
small scale spectrum with 77 ft 3.5, which dominates the cal¬ 
culation of rms curvature, should be visible in the scaling re¬ 
lation for Cims if the resolution is better than A^min = 0.01 
mm. Our resolution is almost always better than this. Because 
we only see the steeper spectrum in the Cnns scaling (Eig. 5), 
the Kelvin spectrum near the resolution might be steeper than 
77 = 3 . 5 , or alternative because the assumption |w'(z)| <C 1 is 
likely not valid, the rms curvature grows slower than predicted 
by Eqs. (4). Also, the above approximation is rather sensitive 
to the value of <I>. By using <I> = 0.05 the critical resolution 
increases to A^niin = 0.0006 mm, which is the best resolution 
used. Now, the cusp contribution should always dominate. 


V. DISCUSSION 

In more complicated tangles, where vortices experience a 
large number of reconnections, the rms curvature can remain 
high only if the average time between successive reconnec¬ 
tions is smaller than the time scale for mutual friction damp¬ 
ing. This implies that the temperature must be low enough, 
as observed in simulations depicting decaying counterflow 
turbulence^. By considering only a single reconnection event 
we are still far from explaining the dissipation in more com¬ 
plicated tangles where reconnections occur continually. The 
future question might be to explain what happens to existing 
Kelvin waves when a new reconnection occurs. Strong local 
stretching caused by a reconnection is likely to damp the pre¬ 
existing Kelvin waves near the reconnection site. 

Our method of determining the Kelvin spectrum from the 
decay of Crms, i e., obtaining 77 from the fitted value of the de¬ 
cay exponent j3, could also be applied to more complicated 
vortex tangles. If one needs to avoid reconnections one could 
apply the method separately for each vortex. One may use 
rather high values of mutual friction and simulate only over a 
short time window to obtain the scaling law. Actually, with 
large enough mutual friction the effects from the cascade are 



FIG. 8 . (Color online) Decay of the rms curvature in the case of 
a straight vortex when a = 0.1 and Lr = 1 mm. The initial Kelvin 
spectrum, shown in the inset, is the Vinen spectrum with 7] = 3 and 
A = O.OIL; at low k’s, while at high fc’s the spectrum is parametrized 
by 77 = 4 and A = O.lLj. The early time behavior (blue curve) is 
determined using 65536 points on the vortex, while 8192 points is 
used to obtain the late time asymptotic (red curve). The dashed lines 
are the analytical results, Eq. (12), corresponding to both spectra. 


minimized and one obtains more accurately the Kelvin spec¬ 
trum which describes the initial state of the vortex. However, 
with large a we cannot immediately tell whether the spectrum 
is due to a cascade or due to other means like excitations from 
reconnections. If the cascade should become important at low 
enough temperatures, it should then result in the decay law 
of Eq. (12) not being satisfied anymore. This change can be 
verified by repeating the decay analysis at these lower temper¬ 
atures. 

A nonconstant decay exponent j3 is likely to originate from 
a Kelvin spectrum where 77 is not universal. This is supported 
by the simulations, similar to the ones presented in Sec. Ill A, 
where we have occupied a straight vortex by using the Vinen 
spectrum at low k’s (m < 100) with amplitude A/L, = 0.01, 
while using a steeper spectrum with 77 = 4 at high k’s with 
A/L; = 0.1, such that the spectrum remains continuous. The 
resulting decay curve has a wide crossover region where p 
changes but both asymptotics follow nicely the theory predic¬ 
tion, Eq. (12). This is illustrated in Eig. 8 . 

The initial approximately constant value of the rms curva¬ 
ture, seen, e.g., at small values of at in Eig. 1, illustrates that 
the numerical resolution determines when the simulation is 
numerically in the zero temperature limit. In this limit even 
the smallest resolvable scales are not smoothened out by mu¬ 
tual friction. Erom the decay of the smallest scale Kelvin 
waves one may derive that if a At (A^ies)^/^ then the ef¬ 
fects from mutual friction are negligible within times of order 
A7. This is consistent with the observations by Kondaurova et 
al.^ where their simulations at T < 0.5 K are essentially the 
same as those conducted at T = 0 . 

The presence of the Kelvin-wave cascade in simulations has 
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been a controversial topic. The concern concentrates around 
the computational resolution required to resolve the excita¬ 
tions at short length scales in the zero-temperature limit. In 
Ref. [15] Kivotides et al. study the reconnection of four rings 
in r = 0 simulations and report excitation of Kelvin waves on 
different scales, including the appearance of a cascade. The 
later conclusion is not in agreement with our findings. How¬ 
ever, in the simulations of Ref. [15] the extremely crinkled 
shape of the vortex after the reconnection differs from simu¬ 
lations conducted later^’*® where more effort has been spent 
on the numerical stability and conservation of energy. See 
also Ref. [9] for numerical challenges that appear when Kelvin 
waves are being identified. 


VI. CONCLUSIONS 

Based on vortex filament simulations for a single recon¬ 
nection event in superfluid helium, we have developed an ex¬ 
planation for the scaling relation seen for the rms curvature 
as a function of the numerical resolution. In contrast to sug¬ 
gestions by Kondaurova et al.^, our model does not involve 
the Kelvin cascade but originates from the excitations cre¬ 
ated by a reconnection cusp. The distribution of Kelvin waves 
with different wavelengths, which becomes visible after the 


cusp relaxes, can be understood to arise from the prereconnec¬ 
tion dynamics where the minimum distance between vortices 
necessarily sweeps all length scales down to the core scale. 
Therefore, the cusp is built from Kelvin waves of different 
scale, limited only by the numerical resolution, that redis¬ 
tribute when the cusp relaxes. 

Thus we conclude that the response to a single reconnection 
can take place without the excitation of a Kelvin-wave cas¬ 
cade. This statement was explicitly tested in vortex filament 
calculations for mutual friction dissipation down to a > 10^^ 
and additionally at T = 0. Instead the reconnection cusp leads 
to an exponential decay of the calculated rms curvature Cms 
which can be associated with a Kelvin spectrum n{k) 
in the range rj ~ 3...3.5. We believe that this indicates that 
a universal Kelvin wave spectrum is not necessarily to be ex¬ 
pected as a response to a single reconnection event. 
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